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Abstract. In this paper, we test the spherical collapse of a unified dark fluid (UDF) which has constant 
adiabatic sound speed. By choosing the different values of model parameters Bs and a, we show the 
nonlinear collapse for UDF and baryons which are considered for their formation of the large scale structure 
of our Universe. The analyzed results show that larger values of a and Bs make the structure formation 
faster and earlier. 
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1 Introduction 

A unified dark fluid (UDF) model [Tl[^|^|il[Hl[^[71[Rl|^[Tni 
fTTlfT2] . which unifies cold dark matter and dark energy 
as a unique component, was considered as an alternative 
model to explain the late time accelerated expansion of our 
Universe [HlfM] . Usually the UDF behaves like the cold 
dark matter at the early epoch and like the dark energy at 
the late time. On the background level, the UDF model 
can mimic the ylCDM model very well. And the cosmic 
observations from the type la supernova (SN) , the Baryon 
acoustic oscillation (BAO), and the full information of the 
cosmic background radiation (CMB) can not distinguish 
the ylCDM model from the UDF model, for example please 
see [8l[8|l9]. 

But can the UDF form the observed large scale struc- 
ture of our Universe? The primordial quantum perturba- 
tions seed the large scale structure of our Universe. And 
the small perturbations grow into the large scale struc- 
ture via gravitational force due to the local over density 
contrast and gravitational instability. So for any realis- 
tic cosmological model, it is demanded to describe the 
whole background evolution history and the formation of 
the observed large scale structure. It is also the situation 
for the UDF models without exception. If a cosmological 
model cannot describe the observed large scale structure, 
it would be ruled out through a comparison of the the- 
oretical calculation and the cosmic observations. When 
one investigates the evolutions of density perturbations 
of UDF, the non-linear stages of perturbations are not 
avoided. However, the study on the non-linear evolutions 
is expensive when one uses the hydrodynamical/A^-body 



numerical simulation. Before the expensive investigation, 
we should test the viability for the formation of the large 
scale structure via a simple model or framework. Fortu- 
nately one can investigate the non-linear stages of collapse 
in the framework of spherical " top-hat" collapse and leave 
the cumbersome hydrodynamical/iV-body numerical sim- 
ulation alone. 

Actually, in the last few years, the generalized Chaply- 
gin gas (gCg) model [ll[51ini[ZllH] , which is an extension of 
the standard Chaplygin gas model [TS] , as a famous UDF 
model was discussed extensively. It covers the range from 
the background evolution and the linear perturbations to 
the non-linear collapse [TB1IT7] in the framework of spher- 
ical collapse model. However when one studies the non- 
linear perturbations, one has to overcome the so-called 
averaging problem [51118). The problem comes from the 
fact that (p) ^ p{{p)) for the gCg model [1] 
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where A and P are positive model parameters, one can 
check that 



{p)^-{A/p^)^-A/{pf=p{{p)), 



(2) 
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in the case of /? / 0. However, it is not the issue for a 
model which has a linear relation between p and p, i.e. p = 
ap — A. This model has a constant adiabatic sound speed 
and has been investigated extensively in the last few years 
2,3,fTU]. It would be better to investigate the non-linear 
perturbation evolution in this model due to the escape 
from the averaging problem. But whether can this UDF 
form the large structure of our Universe? And what are the 
non-linear perturbation evolutions for this UDF? What 
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are the elFects of the model parameter on the collapse? 
To answer these questions is the main motivation of this 
paper. But instead of using the hydrodynamical/A'^-body 
numerical simulation, we will check the possibility of the 
collapse or the large structure formation in the framework 
of spherical collapse for the UDF with a constant adiabatic 
sound speed. 

The paper is structured as follows. In section [2j we 
give a brief review of the UDF with a constant adiabatic 
sound speed. We present the basic equations for spherical 
top-hat collapse of fluid in section [3] The main results are 
summarized in section H) Section [5] is the conclusion. 



2 A Brief Review of Unified Dark Fluid with 
Constant Adiabatic Sound Speed 

In this section, we will give a brief review of a UDF model 
which has a constant adiabatic sound speed (CASS) cl — 
a, for the details please see the Ref. From the definition 
of the adiabatic sound speed 



dpd 
dpd 
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dpd 
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dwd 
dpd 



Wd = a, 



after integration, one obtains 



Wd 



A 



(3) 



(4) 



where A and Wd are the integration constant and the equa- 
tion of state (EoS) of the UDF respectively. In terms of 
a, one has the pressure of the UDF pd — apd — A. The 
energy density and the EoS of this UDF can also be recast 
into the following forms 

= P<io{(l-Bs) + S,a-3(i+")}, (5) 
(l + a)(l-B,) 



Wd = a 



(l-S,)+B,a-3(i+")^ 



(6) 



where in the range < < 1 and a in the range 
< a < 1 are the model parameters. When Bg is zero, 
this UFD is a cosmological constant. If a = is respected, 
this UDF is the standard ylCDM model at the background 
level. And when Bs = 1 and a = are set simultaneously, 
the cold dark matter is arrived. In a spatially flat FRW 
Universe, the Friedmann equation is given as 



i?^ = i?o |(i - - a,: 



(l-i3,,)+S.a"'^'+"^ 



(7) 
}, (8) 



where H is the Hubble parameter with its current value 
Ho = 100/ikm s"^Mpc"\ and f2i [i = b,r) are dimension- 
less energy density parameters for baryon and radiation 
respectively. 

In Ref. U] , this model was tested by using SN la, BAO 
and the full information from WMAP-7yr data sets, where 
the model parameter space was obtained. The result shows 



that at the background level the currently available cos- 
mic observations cannot distinguish the UDF model from 
the ylCDM model. It means that they have almost the 
same background evolution history. However the dynamic 
evolutions would be strongly different even they have the 
same background evolution. One should use the large scale 
structure information to test the viability of UDF models 
as additions. Before doing that, at first we should check 
the possible collapse under this model. In the following 
section, we will use the spherical collapse model to inves- 
tigate the non-linear evolution of the UDF perturbations. 



3 Basic Equations for Spherical Top-hat 
Collapse of Fluid 

The spherical collapse (SC) describes the evolution of a 
spherically symmetric perturbation embedded in a ho- 
mogenous static, expanding or collapsing background. In 
this paper, for simplicity, we consider the spherical top- 
hat collapse (STHC) model which is characterized by the 
perturbed region with constant density [TB] . With this as- 
sumption in hand, the perturbation only depends on the 
time. And one doesn't need to consider the gradients in- 
side the perturbation region. 

In the STHC model, the basic equations which depend 
on local quantities can be written as 



Pc = -3h{pc+Pc), 
r AttG 



'^{Pc, +3pcJ, 



(9) 

(10) 



where h = r/r and r are the local expansion rate and the 
local scale factor respectively. Here the perturbed quanti- 
ties Pc and Pc are defined was 



Pc = p + Sp, 
Pc=p + Sp, 



(11) 
(12) 



where p and p are background quantities which respect to 
the background evolution equations 

(13) 

(p, + 3p,), (14) 



p=-3Hip + p) 
a _ 47rG >^ 
a 



where H = a/a is the Hubble parameter which relates to 
local expansion rate in the STHC model by [T7] 



h = H 



3a ' 



(15) 



where 9 is the divergence of the peculiar velocity i.e. 
6* EE V • ^. 

The governing equations which describe the dynamical 
evolution of density contrast 5i — {Sp/ p)i and 9 are given 
in the following form P^ITT] 

S, = -SHici -w,)S,-[l + w, + {l + cl)6,]-, (16) 

a 
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where Cg. = {5p/5p)i is the square of the effective sound 
speed of energy component i. One can recast the Eq. (jl6p 
and Eq. (|17p into the equations which evolve with respect 
to the scale factor a 



-5(4 

a 



Wi)5i ~ [1 

'2 m 



'Sa^H 



+ (1 ■ 



34), 



,(18) 
(19) 



where we have used the definition f2i ~ SirGpi/SH^ . 

Observing the above dynamical evolution equations, 
one can find that the important involved quantities are 
the EoS Wc of the collapsing region and the effective sound 
speed Cg. One can read off the EoS Wc from its definition 

m 



P + 6p 



w 



+ c 



(20) 



p + Sp 1 + S "l + S 

So, the most important quantity is the effective sound 
speed Cg. For the CASS model, the effective sound speed 
is given as 

2 Sp Pc~P 
dp Pc- P 
By using the relation p — ap ^ A, one has 



(21) 



japc) -A-{ap-A) 
Pc- P 



(22) 



without surprise. 



Model 


a 


Zta 


Sb{zta)/5d{zta) 


a 





0.0678 


1.240 


b 


10"^ 


0.111 


1.211 


c 


10'^ 


0.138 


0.689 


d 


10-1 


0.940 


0.010 



Table 1. Models for the STHC model, where the values of a 
are small positive values due to the constraint from background 
evolution history. The redshift zta denotes the turnaround red- 
shift when the perturbed region begin to collapse. 



10'°- 
10* 
10"* 

10" 

100 

1 







10» 
10'* 
10" 
100 

1 





4 Perturbation Evolutions of Density Contrast 

Now, one can investigate the evolutions of density con- 
trast of different energy component numerically in a spa- 
tially flat FRW Universe, where the relevant cosmological 
model parameters are borrowed from the results obtained 
in Ref. 3: Ho = 71.341km s"^Mpc"\ il^ = 0.956 and 
f2b — 0.044. Here wc just consider the baryon and UDF 
due to the possible formation of the large scale structure. 
We solve the differential equations by using the software 
Mathematica and setting the initial conditions (ICs) Sd 
and 5b at the redshift z = 1000 which is the same redshift 
used in Ref. [T6l. 

To investigate the effects of the model parameter a to 
the spherical collapse of baryon and UDF, we fix the other 
relevant model parameter Bg to its central value = 
0.229 and initial conditions 6d{z = 1000) = 3.5 x 10-^ and 
5b{z = 1000) = 10^^ by varying the model parameter a in 
the range from to 1. We show the calculated results in 
Table[TJ where the redshift zta is defined as the turnaround 
redshift when the perturbed region begin to collapse [TBI. 
The evolutions of the density perturbations for baryons 
and UDF are shown in Figure [T] The vertical parts of the 
curved lines denote the collapse of the perturbed regions. 
The corresponding turnaround redshifts can be found in 
the Table [TJ Higher values of a result in values of Wc closer 
and higher to during the collapse as shown in Figure 
[21 Then the perturbations collapse earlier for the larger 
values of a as shown in Figure [TJ 



0.01 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

z 

Fig. 1. The evolutions of density perturbations with respect 
to the redshift for different models [T] The top and bottom 
panels are for baryons and UDF respectively. Where the or- 
ange, dashed, dotted and green curved lines are for a — 
0,10-^,10-2,10-1 respectively. The horizon line denotes the 
limit of linear perturbation, i.e. 5 — 1. The vertical parts of 
the curved lines denote the collapse of the perturbed regions. 
The corresponding turnaround redshifts can be found in the 
Table [l] 



It's time to show the influence of ICs on the evolution 
of the perturbations. Here we adopted different values of 
ICs for the UDF as shown in Table [H The corresponding 
evolutions of density perturbations are shown in Figure 
[3] One can see that smaller value of initial conditions can 
result in later collapse as expected. 



In the above, we have fixed the values of model param- 
eter Bs to its central value obtained by cosmic background 
evolution constraint. We should also investigate the effects 
to the collapse of the UDF of this model parameter. As 
pointed in our previous work |3] , the model parameter B^ 
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-0.25 - 



-0.30 ' 1 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

z 

Fig. 2. The evolutions of Wc and Wd with respect to the redshift 
z for different models [T] The top and bottom panels are for Wc 
and Wd respectively. Where the orange, dashed, dotted and 
green curved lines are for a = 0, 10"^, 10"", 10"^ respectively. 
Larger values of a result in values of Wc closer and larger to 
during the collapse. 



Model 


a 


ICs (xlO-^) 




Sb{zta)/Sd{zta) 


a 





3.500 


0.0678 


1.240 


b 


10"^ 


3.490 


0.111 


1.211 


c 


10"" 


3.460 


0.122 


0.862 


d 


10-1 


3.210 


0.778 


0.00847 



Table 2. Models for the STHC model, where different initial 
values of UDF are adopted. The corresponding turnaround red- 
shifts are summarized. 



effectively takes the role as the ratio of cold dark matter 
in our Universe. Then it is easy to understand that the 
larger values of Bs will increase the turnaround redshift 
Zta- To see that, we show the zta with respect to Bg for 
different values of a in Figure SI 



Through the above calculation and analysis, one can 
clearly see the possible large scale structure formation in 
this CASS model and understand the effects of the model 
parameters to the evolutions of the density perturbations. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Fig. 3. The evolutions of density perturbations with respect 
to the redshift for different models [S] The top and bottom 
panels are for baryons and UDF respectively. Where the or- 
ange, dashed, dotted and green curved lines are for the models 
a = 0, 10"^, IQ-", IQ-i respectively. The horizon line denotes 
the limit of linear perturbation, i.e. 5—1. The vertical parts of 
the curved lines denote the collapse of the perturbed regions. 
The corresponding turnaround redshifts can be found in the 
Table [1 



0.18 r 
0.17 : 
0.16 : 




0.12 : - 
0.11 : 

0.10 E — ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' — : 

0.20 0.22 0.24 0.26 0.28 0.30 

Fig. 4. The turnaround redshift Zta with respect to model 
parameter a for different values of a. Where the orange solid 
line is for a = 0.01 and the dashed line is for a = 0.009. 
For plotting this figure, the initial conditions Sd{z = 1000) — 
3.5 X 10"^ and 6b{z = 1000) = 10"^ are adopted. 



5 Conclusion 

In this paper, we checked the possibility of the large scale 
structure formation for a UDF model which has a con- 
stant adiabatic sound speed (CASS) in the framework 
of spherical collapse. By choosing the different values of 
model parameters a and B^, we show their effects on the 



non-linear perturbation evolutions. The results show that 
larger values of a and Bs can make the density perturba- 
tions collapse earlier. After the investigation, we can con- 
clude that it is possible to form large scale structure in 
the CASS UDF model via the spherical collapse. And the 
next step, we should study the hydrodynamical/iV-body 
numerical simulation and compare the simulated results 
with the observed large scale strutter of our Universe. 
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